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ABSTRACT 

We discuss general relativistic effects in the steady-state neutrino-driven 
“wind” which may arise from nascent neutron stars. In particular, we generalize 
previous analytic estimates of the entropy per baryon S, the mass outflow rate 
M, and the dynamical expansion time scale Td yn - We show that S increases 
and Tdyn decreases with increasing values of the mass-to-radius ratio describing 
the supernova core. Both of these trends indicate that a more compact core 
will lead to a higher number of neutrons per iron peak seed nucleus. Such 
an enhancement in the neutron-to-seed ratio may be required for successful 
r-process nucleosynthesis in neutrino-heated supernova ejecta. 


Subject headings: nuclear reactions, nucleosynthesis, abundances — supernovae: 
general — equation of state — relativity 

The production site of the r-process elements (Burbidge et al. 1957; Cameron 1957) is 
a longstanding problem (Mathews & Cowan 1990). One of the most promising candidate 
sites for r-process nucleosynthesis is the neutrino-heated ejecta from the post-core-bounce 
environment of a Type II or Type lb supernova (Meyer et al. 1992; Woosley & Hoffman 
1992; Takahashi, Witti, & Janka 1994; Woosley et al. 1994). These r-process calculations, 
though promising, cannot reproduce a solar system r-process abundance pattern without an 
artificial increase in the neutron-to-seed nucleus ratio over that predicted in hydrodynamical 
calculations and simple wind models (see, e.g., Hoffman, Qian, & Woosley 1996; Meyer, 
Brown, & Luo 1996). 

Qian & Woosley (1996, hereafter QW) used a simple model of the neutrino-driven 
wind to obtain both analytic and numerical estimates of quantities upon which the 
nucleosynthesis abundance yield depends. These quantities include the electron fraction Y e , 
the entropy per baryon S, the mass outflow rate M, and the dynamic expansion time scale 
r dyn —all of which are important in setting the neutron-to-seed nucleus ratio prior to the 
epoch of rapid neutron capture. 


2 


In this Letter we generalize the analytic derivations in QW to include general relativistic 
effects. Such effects are of potential interest in light of “soft” nuclear matter equations 
of state—involving, for example, kaon condensation (Thorsson, Prakash, & Lattimer 
1994)—which could lead to very compact supernova cores and even black holes as supernova 
remnants (Bethe & Brown 1995; Woosley & Timrnes 1996). While QW reported two 
sample numerical calculations involving post-Newtonian corrections, they did not present 
corresponding analytic calculations. Numerical calculations reported in QW suggest that 
general relativistic effects go in the direction of making conditions more favorable for the 
r-process. We here allow for general relativistic effects in semi-analytical calculations, 
including effects not included in the numerical calculations by QW, namely the redshift of 
neutrino energies and bending of the neutrino trajectories. 

In the analytic approximations performed by QW, the calculation of the quantities S, 
M , and rd yn essentially decouples from the calculation of Y e . We will present semi-analytic 
estimates for S, M, and Td yn ; general relativistic effects on Y e have been considered 
in another paper (Fuller & Qian 1996). Unless units are explicitly given, we take 
h = c = k = G = 1. 


The wind equations in Duncan, Shapiro, and Wasserman (1986) and QW can be 
generalized to allow for relativistic outflow velocities and general relativistic effects in a 
static Schwarzschild spacetime. A detailed derivation and discussion will appear elsewhere 
(Cardall & Fuller 1997, hereafter CF). We here simply present the general relativistic 
analogs of Eqs. (24-26) of QW, which give the radial evolution of the velocity, “flow energy” 
per baryon £fl ow , and entopy per baryon: 
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In these equations, v is the outflow velocity as measured by an observer at rest in the 
Schwarzschild spacetime; v s = ( TS/Sm^) 1 / 2 is the sound speed; T is the temperature; S is 
the entropy per baryon; mjv is the baryon rest mass; M is the mass of the supernova core; 
q is the heating rate per unit mass; and y = [(1 — 2M/r)/(l — u 2 )] 1 ^ 2 . These equations 
assume radiation dominated conditions, i.e. that the pressure and energy density are 
dominated by relativistic particles; a steady state outflow; and that the supernova core is 
the dominant source of gravity. 
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The specific net heating rate q includes a number of processes. Heating processes 
include (anti)neutrino absorption on nucleons, ve scattering, and vv annihilation. The 
major cooling processes are e ± capture on free nucleons and e ± annihilation. All but the 
last of these processes depend on neutrino luminosities, energy moments of the neutrino 
distribution functions, and geometric factors involving the maximum neutrino deviation 
angle from the radial direction at a given radius. Expressions for these rates are given in 
QW. General relativistic effects can be introduced through redshift factors in the energy 
moments, and redshift and time dilation factors in luminosities ([Fuller fc Qian 1996|) ; and 
by suitably altering the geometric factors to account for curved neutrino trajectories (|CF| ) . 


Before proceeding to semi-analytic estimates, we briefly discuss initial and boundary 
conditions. Treatment begins at the surface of the supernova core, taken to be Schwarzschild 
coordinate radius R. As in QW, the initial entropy is taken to be Si = 4 (S > 4 assures 
radiation dominated conditions). The initial temperature is obtained by equating the 
heating and cooling rates at R. In QW neutrino absorption could be taken as the sole 
heating process. This allowed a fairly simple expression for to be obtained in terms of 
the v e luminosity, the supernova core radius R , and the ratio of the second and first energy 
moments of the v e spectrum (e 9f , = (E? e )/(E Pe )). For more compact cores, however, the 
other heating rates become more important at the neutron star surface and it becomes 
impossible to obtain an analytic expression for 7) in terms of the above quantities. Therefore 
we choose the v e luminosity to be 10 51 erg s -1 , M = IAMq, and e Pe = 20 MeV, and 
numerically solve for 7j for various values of 2 M/R. 


Neutrino heating becomes negligible by T 0.5 MeV, where efl ow and S reach their 
final values, after which they remain constant. These final values will be denoted by a 
subscript /. We assume that a boundary is provided by the supernova shock at some large 
radius, where the temperature is X),, and where the flow asymptotes to small velocity. We 
then have from Eq. ([|) 



between T — 0.5 MeV and T b . 


We now obtain an approximate expression for Sf. From Eq. (^;), 
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where T eff and S e s will be defined below. From Eq. (|3|), 

Sf ~ S f -Si = m N f f ^ — . (7) 

Jr T vy 

Combining these last two equations, taking In (1 + T b Sf/rriN ) ~ T b Sf/m,N, and using the 
definitions (QW) S eS ~ Sf/ 2 and T eff = [J (vy)~ 1 q dr] / [ /{vyT)~ l qdr\ , we obtain 

„ mjv In (1 — 2M/R) 

S f ~-— 77 :-- — 1 

J 2T efr 

The quantity T eff is an average temperature, weighted by the heating rate. In QW it 
is approximated as the temperature at which the heating rate is a maximum. Using 
only neutrino absorption and annihilation, they obtain T e g = 6 -1 / 6 T t . As mentioned 
previously in connection with the determination of the initial temperature T l} addition 
of the other heating and cooling processes makes analytic progess impossible. To obtain 
approximate numerical results we will employ the expression T e g = 6 -1 / 6 R as obtained by 
QW. Figure 1 is a plot of S as a function of 2 M/R, obtained from Eq. (H) with T b = 0.1 
MeV. The circle shows the numerical result obtained in QW for their model 10B with 
post-Newtonian corrections. The agreement is excellent. 
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Next we estimate M. Baryon mass conservation in the Schwarzschild geometry gives 
M = 47rr 2 p b vy = constant, where p b is the baryon mass density &■ As in QW, we 
evaluate this at r e fj « R to obtain M = Ai tR 2 p bje s ( vy) e s . The quantity p bfi s can be derived 
from the entropy per baryon SOr = (ll7r 2 /45)(m 7 vT e 3 ff /pb ie fj). To estimate (vy) e g, we use Eq. 
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which can be solved for (vy) e g. The quantity g eff is obtained by evaluating the net heating 
rate q at temperature T e and radius r e s ~ R. We follow QW and estimate the temperature 
gradient by assuming approximate hydrostatic equilibrium near radius R ; the result is 
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Now M can be calculated using the previously obtained values of S. A plot of M as a 
function of 2 M/R is given in Figure 2. (Note the kink at R = 3M; below this radius, there 
is a maximum angle of deviation from the radial direction beyond which massless particles 
cannot escape to infinity.) 


The dynamic expansion time scale 7d yn is defined to be 
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where the subscript / means that the quantity is evaluated at T — 0.5 MeV. The dynamic 
expansion time scale is closely related to the proper time a fluid element experiences in 
going from Tj = 0.5 MeV to, say T 2 = 0.2 MeV: 
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In this case, we estimate the temperature scale height by approximating the boundary 
condition in Eq. (|4|) as TSf/rriN ~ M/r. Then 
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The quantity (vy)f can be obtained from the previously computed quantities M and Sf. A 
plot of r dyn as a function of 2 M/R is given in Figure 3. 


Neutron-to-seed ratios on the order of 100 or greater are necessary for a successful 
r-process, and this can be achieved under a variety of combinations of S, Td yn , and Y e 
(Hoffman et al. 1996; Meyer et al. 1996; Meyer & Brown 1997). We have here confirmed 
that general relativistic effects increase S and reduce Td yn , as seen in selected numerical 
calculations in QW. Both of these trends lead to a higher neutron-to-seed ratio. On the 
other hand, Fuller & Qian (1996) have shown that general relativistic effects tend to 
increase Y e because of the differential redshift of the u e and v e emitted from the supernova 
core. This differential redshift decreases the difference in the average energies of the v e and 
v e populations, driving Y e larger and so closer to 0.5. This is probably unfavorable for a 
high neutron-seed-ratio, since the general trend is that a higher Y e requires a higher S to 
obtain a given neutron-to-seed ratio. The magnitude of the differential redshift effect on Y e 
is uncertain because of its dependence on unknown details of the nuclear equation of state. 


On the whole, however, general relativistic effects are likely to increase the neutron-to- 
seed ratio. For one thing, an increased Y e due to differential neutrino redshifts is not always 
bad: examination of Figure 10 and Table 5 of Hoffman et al. (1996) shows that the entropy 
requirements to obtain a given neutron-to-seed ratio actually become less severe as Y e gets 
very close to 0.5. More important than this, however, is the fact that salutory effects on 
S and rd y n can more than compensate for deleterious differential redshift effects on Y e . 

For example, consider the case of 2 M/R = 0.68, for which our semianalytic estimates give 
Td yn ~ 0.004 s and S ~ 255. In Table 5 of Hoffman et al. (1996), for “expansion time” 
0.005 s (which corresponds to Td yn — 0.004 s), an entropy of only 147 is needed to obtain a 
neutron-to-seed ratio ^ 100, even for the “worst case” value of Y e . 

It should be pointed out that the case we just considered, 2 M/R = 0.68, (just) 
violates the causality limit R > 1.52 x 2 M imposed on any equation of state allowing a 
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stable neutron star. However, we only consider this particular case for comparison with 
Hoffman et al. (1997). Unfortunately, their Table 5 does not contain values of expansion 
time between 0.005 s and 0.025 s. We fully expect that there are values of expansion time 
between 0.005 s and 0.025 s for which 2 M/R does not violate the causality limit, and for 
which S is sufficiently high to provide a neutron-to-seed ratio in excess of 100 even for the 
“worst case” Y e . Also, even beyond the causality limit, our calculations may give an idea of 
what happens near an unstable supernova core that collapses to a black hole at relatively 
late times. Of course, the time scale of collapse may be fast enough to grossly violate our 
assumption of static conditions. Reliable exploration of this case would require detailed 
numerical modeling and good knowledge of the equation of state of nuclear matter. 

There are other questions regarding the viability of a highly relativistic supernova core 
as a suitable environment for r-process nucleosynthesis. While a small dynamic expansion 
time scale is conducive to a large neutron-to-seed ratio, it may not allow enough neutron 
capture time to build up the r-process elements from the seed nuclei. Some mechanism 
would have to exist to slow down the flow during the neutron capture epoch (see e.g. 
McLaughlin & Fuller 1997); perhaps the shock used in the boundary condition could serve 
this purpose. Also, we have seen that a more compact core leads to a smaller mass outflow 
rate. One of the attractions of neutrino-heated supernovae ejecta as an r-process site is 
that the total mass loss, together with the estimated Galactic supernova rate, roughly fits 
the observed amount of Galactic r-process material (Meyer et al. 1992). One would not 
want M to become so small that this agreement is ruined. However, even for the extreme 
case of 2 M/R = 0.68 considered above, our estimated M ~ 1.3 x 10 _6 Mqs _1 may be large 
enough to be viable. 


We wish to thank Y.-Z. Qian and J. R. Wilson for helpful discussions. This work was 
supported by grants NSF PHY95-03384 and NASA NAG5-3062 at UCSD. 
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Fig. 3.— The dynamic expansion time scale as a function of the supernova core Schwarzschild 
radius divided by the core radius. The circle is from the Qian & Woosley (1996) numerical 
calculation of model 10B with post-Newtonian corrections. 



